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Spectrum of an open disordered quasi-two-dimensional electron system: 
strong orbital effect of the weak in-plane magnetic field 
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The effect of an in-pfane magnetic fieid upon open quasi-two-dimensionai eiectron and hoie systems 
is investigated in terms of the carrier ground-state spectrum. The magnetic fieid, classified as weak 
from the viewpoint of correlation between size parameters of classical electron motion and the gate 
potential spatial profile is shown to efficiently cut off extended modes from the spectrum and to 
change singularly the mode density of states (MDOS). The reduction in the number of current- 
carrying modes, right up to zero in magnetic fields of moderate strength, can be viewed as the cause 
of magnetic-field-driven metal-to-insulator transition widely observed in two-dimensional systems. 
Both the mode number reduction and the MDOS singularity appear to be most pronounced in the 
mode states dephasing associated with their scattering by quenched- disorder potential. This sort 
of dephasing is proven to dominate the dephasing which involves solely the magnetic field whatever 
level of the disorder. 

PACS numbers: 71.30.+h, 72.10.(1, 72.15.Rn, 73.23.b, 73.50.-h, 73.40.Qv 
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I. INTRODUCTION 

The apparent metallic state widely observed in two- 
dimensional (2D) systems of Si MOSFET type as well 
as in GaAs/AlGaAs heterostructuresi^ obviously con- 
tradicts the well-known one-parameter scaling theory 3 
and as yet has not received generally accepted theoreti- 
cal explanation. The existence of such a state is mostly 
believed to result from Coulomb interaction of carriers, 
which is rather strong in the systems of low electron and 
hole density. Estimations of this interaction indeed can 
cause the surmise that just this interaction should lead to 
quite strong dephasing effect upon electrons which oth- 
erwise would be localized due to scattering by the disor- 
der potential, thus residing in coherent states. However, 
at the present time the lack of a comprehensive theory 
for Coulomb interaction in solids precludes from making 
certain conclusions about its predominant role in form- 
ing the metallic ground state of 2D electron systems. In 
particular, a considerable challenge in this connection is 
presented by commonplace observations of the dephasing 
time saturation in different systems, including 2D ones, 
at temperatures very close to zero (see, e. g., Ref. and 
numerous references therein). 

Besides the unexpected conduction state of quench- 
disordered 2D systems, the metal-insulator transition 
(MIT) is normally observed there, which currently has 
not been proven unambiguously to be determined exclu- 
sively by the level of the disorder. No less puzzling is also 
the abnormally large response of 2D electron and hole 
systems to the relatively weak in-plane magnetic fields 
which is known to significantly suppress the metallic be- 
haviour of the carriers and even to drive the system into 
the insulating regime &L&& 

In so far as the electrons confined to move in a nar- 
row near-surface potential well are weakly coupled to the 
in-plane magnetic field through their orbital degree of 



freedom, it is widely believed that such a field promotes 
localization of carriers, and thus the MIT, mainly due to 
strong spin-related effects OL1LH The relatively rare pa- 
pers where orbital coupling was analyzed by taking into 
account the finite width of potential wells forming real 
two-dimensional systems are not rated by now as fully 
convincing. Specifically, the relatively simple model sug- 
gested in Ref. [lj do not exhibit sufficiently abrupt tran- 
sition between metallic and dielectric regimes, whereas in 
Ref. [lj only corrections caused by weak localization of 
carriers are studied, which can hardly serve as the con- 
clusive proof for the physical mechanism of the observed 
effects. 



Previously in Refs. Il5ll6l one-particle theory capable 
of explaining the metallic ground state as well as MIT 
in disordered 2D systems not subjected to magnetic field 
was developed starting from basic positions essentially 
different from those of scaling theories. Specifically, the 
conductance of a strictly 2D— and a quasi-2Di& system 
was calculated in terms of quantum states pertinent to 
a perfect finite-size open system of waveguide geometry. 
In this approach, the metallic value of the conductance 
is bound up with the primary existence of coherent ex- 
tended waveguide modes rather than one-particle elec- 
tronic states originally localized by the disorder. Energy 
levels of these collective mode states can be widened by 
quenched disorder provided that scattering is ensured be- 
tween extended modes having different longitudinal en- 
ergies. This type of scattering can be viewed (math- 
ematically) as inelastic, although it is physically pro- 
vided by a static random potential. In the suggested 
approach, all extended modes other than the particular 
one, if any, can be regarded as the dephasing bath. With 
gradual strengthening of the disorder, the conductance 
transforms from its ballistic value in a perfect system, 
which equals the number of extended modes times the 
conductance quantum, to the diffusive value coincident 



with standard Drude conductance if the system possesses 
the number of extended modes noticeably greater than 
unity. 

In th e pres ent work, we apply the mode approach 
of Refs. 1151161 to examine the influence of the in-plane 
magnetic field upon spectrum of the electrons restricted 
to move in a planar, yet three-dimensional, open quan- 
tum well. It will be shown that at large values of the 
transverse aspect ratio of such an electron waveguide 
even a rather weak magnetic field can significantly af- 
fect the electron spectrum. This appears in considerable 
magnetic-field responsivity of the number of extended 
modes, the latter being normally identified as conducting 
channels, as well as in MDOS sensitivity to the magnetic 
field. The former factor is well known to control the value 
of the ballistic conductance of the confined current carri- 
ers whereas the latter (MDOS) governs substantially the 
mode states dephasing associated with scattering of the 
electrons by the disorder potential. Interestingly, mode 
entanglement solely due to the magnetic field, with no 
disorder whatsoever, leads to electron mass renormaliza- 
tion and does not affect the width of the energy levels of 
the collective electron states. 



II. THE MODEL 

Two-dimensional electron and hole systems in practi- 
cal applications can be modelled, in view of their open 
property in the direction of current flow, as planar quan- 
tum waveguides whose transverse structure is governed 
by lateral (depletion) potential. The exact form of the 
confining potential well is of minor importance for its 
principal application which reduces to the restriction of 
electron transport in the direction normal to the inter- 
facial area and, consequently, to the "transverse" quan- 
tization of the electron spectrum. In this study, we as- 
sume the open quasi-two-dimensional (Q2D) system of 
carriers having the form of three-dimensional "electron 
waveguide" of rectangular cross-section, which occupies 
the coordinate region 



arbitrary except that the thickness H will be assumed to 
serve as a small length parameter which will be specified 
below. The electron system will be thought of as open at 
the ends x = ±L/2 and closed by infinite-wall potentials 
at all lateral boundaries. The magnetic field B is taken to 
point parallel to x-y plane at arbitrary angle with respect 
to x-axes. 

Since main transport coefficients, in particular the con- 
ductance, are expressed in terms of one-particle propa- 
gators of carriers, we will analyze the equation for the 
retarded Green function of Fermi particles with energy 
£p = kp. In the Fermi- liquid approximation it has the 
form 



V-^A(r)) + kp + tO- Tin 



G(r,r') 
= <5(r-r'), (2) 



where $0 = hc/e is the magnetic flux quantum, A(r) is 
the vector potential of the external magnetic field, V(r) 
is the scalar random potential due to, say, impurities or 
rough boundaries of the confining potential well. We 
adopt hereinafter the system of units with H = 2m = 1 , 
to denoting the electron effective mass. 

With the magnetic field gauged so as A(r) = 
(B y z,—B x z,0), the equation J2J) takes the form 



V 2 + k 2 F + iO - V(r) 



Am / d 

$n~ ( ByZ dx~ 



B x z 



d_ 

dy 



$0 



G(r,r')=S(r-r') . (3) 



At this stage it is expedient to go over from the ini- 
tially three-dimensional problem to a set of strictly one- 
dimensional problems, individual for each of the modes. 
Towards this end, as a first step one should carry out 
Fourier transformation of equation ||3J) over the trans- 
verse radius- vector rj_ = {y, z). The appropriate set of 
eigenfunctions has the form 



x€(-L/2,L/2), 
y€[-W/2,W/2] , 
z e [-H/2, H/2] , 



(1) 



as shown in Figure ^ The length L, the width W as well 
as the height H of the model system will be regarded as 




FIG. 1: The geometry of Q2D electron waveguide subject 
to in-plane magnetic field B. 
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where \x = (n, to), with n, to G N, is the vector mode 
index conjugate to the coordinate vector r^. With func- 
tions (0J, equation J3J) is readily transformed to the set 
of coupled equations for mode Fourier components of the 
function G(r,r'), viz. 



d 2 
— + kl + ?,0-V„(x) 



t-^/X/X' \%1 ^ ) 



- ^ tiv,v( x ) G vfj.' {x, x) = (5 M/ y<5(a; - x) . (5) 



Here, 



% = k *-(™Y-(™V 



\W J \ H ) 

is the unperturbed mode energy, 



Mx) = V^(x) 



H 2 
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(6) 



(7) 



is the diagonal-in-mode-indices matrix element of the to- 
tal potential which includes the impurity part V(r) and 
both of the magnetic terms in square brackets of equa- 
tion |J5J, l B = ^/<f>o/27r|B| is the magnetic length. The 
term V^ M (x) in Eq. (JJJ is the diagonal element of the 
mode matrix ||V^i/|| whose components are evaluated as 



Viivix) = dr 
'S 



.(T ±]f t\V{T)\Tr,v) 



(8) 



integration is over cross-section S of the quantum well. 
Off-diagonal mode matrix elements W M „ (x) in Eq. (JSJ also 
include both the disorder and the magnetic-field origi- 
nated potentials, viz. 



ZV(x) =Vy,u{x) + 2iH 
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In Eq. J2J, h is the partial magnetic length given by 
I 2 = $o/2tt|-Bj|, numerical coefficients, which are specifi- 
cal for the geometry of the quantum well, have the form 
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(10a) 



(10b) 



(10c) 



with q = Bi/\Bi\ = sgn Bi. In expressions (|10fl . 
mode indices are designated such that fi = (n, to) and 
f = (ni,mi). 

The potentials V M (x) and U^^x) in Eq. JSJ may be 
thought of as responsible for coherent intra-mode and 
incoherent inter-mode scattering, respectively. We thus 
adopt in this work the approach where interactions of the 
electron system with both the impurities and the mag- 
netic field are exploited on equal footing, that is they 
are treated as the problems of scattering by additive 
static potentials which are basically different in corre- 
lation properties only. 



III. REDUCTION TO ONE-DIMENSIONAL 
DYNAMIC PROBLEMS 

A set of equations JSJ , though describing mode prop- 
agation in one spatial dimension, cannot certainly be re- 
garded as a really one-dimensional dynamical problem by 
virtue of strong correlation of different modes via inter- 
mode potentials © . Formally, this manifests itself in co- 
existence in (jSJ) of purely intra-mode propagators, i. e. 
the Green functions having identical mode indices, and 
inter-mode propagators G vll {x 1 x 1 ) with v ^ [i. 

One can obviate these complications using the method 
suggested in Refs. Il5ll6l In the above papers, non-dia- 
gonal elements of the mode matrix || G M/ v || were proven to 
be expressed, by means of some linear operation, through 
the respective diagonal elements only. Substitution of 
thus represented inter-mode propagators into Eq. JSJ) re- 
sults in the following set of strictly one-dimensional equa- 
tions for intra-mode propagators G MM (x, x'), 



dx 2 

for V/x . 



+ k 2 +iO-V ll (x)-f ll 



G> M (x,x') = 8(x — x') , 

(11) 



Here, 7^ is the operator (integral) potential, well-known 
as the T-matrix in the quantum theory of scattering^i^ 
which acts in the x-coordinate space X. It has the form 



T lx = P l M(t-R)- 1 RP ll 



(12) 



where U and R are the operators acting in the mixed 
mode-coordinate space M M constructed as a direct pro- 
duct of the space X and the truncated mode space which 
incorporates the whole set of mode indices except the 
unique mode index \x. The operator U is specified in M^ 
by its matrix elements 



(x,u\U\x',u') =Uuw(x)S(x-x') 



(13) 



matrix elements of the operator R have the form 

(x,u\R\x',u')=G¥\x,x')U„„,{x') . (14) 

The function G v (x,x') in (jl4|) will be thought of as 
the trial mode Green function which satisfies the equa- 
tion resulting from Eq. ||SJ provided that all inter-mode 
potentials are put identically equal to zero, 



<9x 2 



*s 



*0 - V„(x) 



G^\x,x') = 5(x-x'). (15) 



The operator P^ in lfT2^l is the projection operator whose 
action reduces to assigning the given value fi to the near- 
est mode index of an arbitrary operator standing next to 
it (either to the left or right), without affecting the prod- 
uct in the X space. 

With intra-mode Green functions found from the set of 
equations i|ll|) . all inter-mode propagators are expressed 
via the operator relation 



Gy, 



P v {t - R) 



R-P nG/ifj, 



(16) 



G vfl and G M/X being thought of as matrices in X-space. 
The initially three-dimensional problem thus reduces to 
the set of separate ID equations itTTJl . each representing 
the closed problem provided that trial Green functions 
are independently found from equation (|15|l . 

To analyze the mode states spectrum, i. e. the spec- 
trum of differential operator in Eq. (jl 1|> , it is worthwhile 
to renormalize mode energies in Eqs. I|ll|l and (|15|l by ex- 
tracting from the initial mode energy JfjJ) the non-random 
"magnetic" part of the intra- mode potential V^{x), thus 
defining the new "unperturbed" mode energy 



2 _ ,2 H 

ft fj. l2l i 



B 



1 



<» 



(17) 



In such a way one is led to solve, in place of equations 
(|TT)l and (JTSJ, a couple of different, though equivalent, 
equations, namely 



a 2 

dx 2 



+ ^ +i 0-V^(x)-T, 



G^(x,x') (i8a) 

= 8(x - x') 



and 



— +^ + iO-V^(x) 



G^\x,x')=5{x-x') 



(18b) 



which must be supplied with correct boundary conditions 
at open ends of the system, to be discussed in the next 
section. 



IV. SPECTRUM OF THE MODE STATES 



To examine differential operator QlSaJI . which is of 
principal import for our purpose, one should first solve 
Eq. (|18b|) for the truly one-dimensional trial Green func- 
tion. In the absence of magnetic field, this problem was 
resolved in Ref. [l5| where arbitrary statistical moments 
of the function G^ {x,x') were found in the case of open 
system subject to weak disorder potential by applying the 
averaging technique appropriate for causal-type random 
functionals. The condition for weak impurity scattering 
(WIS), which we assume to hold in this study as well, 
can be cast to the form of the inequality pair 



k' 1 



<^ 



(19) 



were r c is the correlation radius of the random potential, 
£ is the electron mean free path relating to it. 

In the case of non-zero magnetic field the solution to 
Eq. (|18bl) is much more involved than that accomplished 
in Ref. ll5J This will be thoroughly examined in a sepa- 
rate publication, while here we outline the solution along 
with criteria of its applicability. 

Given the magnetic potentials, to adequately take into 
account the open property of one-dimensional system 



governed by equation (|18b|) one should explore this equa- 
tion on the extended x-axis rather than on disordered 
and subject to magnetic field interval C = (—L/2,L/2). 
Being considered on the whole axis, equation (|18b|) 
describes the motion of a quantum particle created 
with energy k 2 at point x' and then propagating 
in two-component scalar potential with bounded sup- 
port. The regular component of this combined po- 
tential is due to the magnetic field. It has the form 
V^\x) = 9(L/2- 1x1)^(1--^), whereas the 

random component, which is of impurity origin, is cov- 
ered by the function V^ an ^(x) = 6{L/2 - \x\)V^(x). 

To perform configurational averaging over the ran- 
dom part of the potential it is worthwhile to express the 
trial Green function in terms of wave functions of causal 
type rather than functions that meet the initially stated 
boundary-value (BV) problem. This may be achieved by 
employing the formula (for the sake of clarity we omit 
mode index fi) 

G(v) ^ x '^ = w ^ + {x) ^- {x ' )e{x ~ x ' ] 

+ ^ + {x')tl}-{x)d{x' -x)] , (20) 

where ip± (x) are two different solutions of homogeneous 
equation 115fl with boundary conditions specified for each 
of them at only one end of the system, viz. x — > ±cx), 
depending on the sign index, W is the Wronskian of those 
solutions. With this representation, the trial propagator 
itself meets, as it must, the initial BV problem. 

The openness of the finite-size system under consid- 
eration implies that far from the source coordinate x' , 
specifically at x — > ±oo, Green function G^ v '(x,x') must 
have the form of outgoing free waves. In view of bound- 
edness of the support of the potentials, at large values 
of |x| functions ip±i x ) have to be taken as 



ift±(x) = c± exp [ ± ik(x — L/2)] 



(21) 



Inside the magnetically biased interval C, in order to 
properly take into account the electron backscattering 
from the potential V^- ran \x), it is worthwhile to seek 
wave functions tp±(x) in the form 

ip±(x) = T:±(x)e l ^ ±x - L ^ - i 1± {x)e~ l ^ ±x - L ^ , 

(22) 

with x 2 specified in (|17|) . Under WIS conditions i|19|l , en- 
velope functions 7r± (x) and 7± (x) in (1221 can be regarded 
as smooth factors in comparison with near-standing fast 
exponentials, which leads to the following coupled dy- 
namic equations, 

± n' ± (x) + ir](x)n ± (x) + (*± (x)l± (x) = , (23a) 
± 7± (s) - iv(.x)l± (s) + C± (a)** (x)=0. (23b) 

Random functions T](x) and C± (x) in Eqs. (|23|l are con- 
structed as normalized packets of spatial harmonics of 



the impurity potential V(x) = Vp^ix), which have the 
form 



T)(X) = 

C±(*) 



1 



x+l 



di 



2xJx-l 2/ 






(24a) 
V(t) exp[2i>c(±x - L/2)] . (24b) 



as weak. The weakness of the impurity scattering implies 
the inequalities (|19fl whereas the magnetic scattering will 
be thought of as weak provided that the requirement is 
met |ft (B) | < 1. From (|25bl) and JT7|) one can make 
sure that in terms of appropriate physical parameters 
the condition for weak magnetic scattering (WMS) may 
be expressed as the inequality 



Spatial averaging in l|24|) is carried out over the interval 
21 of arbitrary length intermediate between small lengths 
>c~ x and r c , on the one hand, and the large scattering 
length (to be determined self-consistently) , on the other. 
In view of these limitations, the "potentials" r\(x) and 
C±{ x ) provide forward and backward scattering of har- 
monics ±>r, respectively. 

By joining the solutions l|22|l and (|21|l at the end points 
of the interval C we arrive at the exact boundary condi- 
tions for the envelopes ir±(x) and 7±(x), viz. 



The quantity 



«1, 



(27) 



7r±(±L/2) =const , 


(25a) 




7 ±(±£/2) =rt B \ ± (±L/2) . 


(25b) 


(V(r)) = 0, 
(V(r)V(r')) = QW(x 



where R c = kpl 2 B is the maximal classical cyclotron ra- 
dius of the electron orbit. It is just from the viewpoint of 
this constraint that we will regard the in-plane magnetic 
field to be weak. 

As far as the impurity scattering is concerned, to do 
the averaging over realizations of the potential V(r) this 
random function will be thought of to possess the follow- 
ing correlation properties, 



x')6(r ± - r' ± ) 



(28a) 
(28b) 



%W = 



(26) 



as it follows from (|22|l . is the amplitude reflection coef- 
ficient from the boundary between magnetically biased 
and unbiased regions. This reflection will be hereinafter 
referred to as "magnetic scattering" associated with the 
above introduced potential V^ reg \x). 

Below in this paper, scattering associated with both of 
the potentials, V^ ran \x) and V^ reB \x), will be regarded 



angular brackets denote configurational averaging. Un- 
der WIS conditions (|19f) . the equalities l)28|) are sufficient 
to adequately accomplish the averaging for rather wide 
class of the random potential statistics, since in this case 
function V(r) may be regarded as approximately Gaus- 
sian distributed^ 

By applying the averaging technique outlined in the 
Appendix the average trial Green function is obtained 
in the following form, which is valid provided WIS and 
WMS conditions hold simultaneously, 



G%\x,x>) 



2x, 



■ exp ■ 



m 2 



,4 y V) 



-(V) 



(M), 



(29) 



Here. 



and 



4 V V) 



4 y V) 



4\ 2 S_ 



W(2x„) 



(30a) 



(30b) 



are the extinction lengths related to forward (/) and 
backward (b) scattering by the potential V lifi (x), \V(q) 
is the Fourier transform of function W(x) from (|28b|l . 
With the result (|29|l , the average square norm of the op- 



erator R from i|12|) can be represented as a sum of "impu- 
rity" and "magnetic" terms, viz. (||R|| 2 ) w (||R(' im P)|| 2 ) + 
(||R( B )|| 2 ), which are estimated as 



<||R«"»>f) ~ JL « 1 , 

<IIR<»¥}~(fT«i. 



(31a) 
(31b) 



Inequalities (|31J) permit simplification of the operator po- 
tential T^ since the inverse operator in (|12fl can be ap- 
proximately replaced with the unit operator. The inter- 
mode potential thus reduces to the relatively simple form, 



t„ « PMG^UP, 



m - 



(32) 



where (/(^ stands for the operator in M^ which is spec- 
ified by the matrix elements of the following form, 



Unlike quasi-local intra-mode potential V, 



hij.\ 



(33) 



the 



operator potential (|3^|l possesses, even in the absence 
of magnetic field, the nonzero mean value. Therefore, to 
apply further a perturbation theory it is reasonable to 
represent this operator as a sum of averaged and fluc- 
tuating parts, i. e. T^ = (7^) + AT^. With regard to 

Eq. |J5J, the mean operator (^) can be splitted (though 
quite conventionally) into the local "impurity" and essen- 
tially non-local "magnetic" terms. The action of short- 
correlated impurity part of this operator reduces to mul- 
tiplication of the mode propagator by the complex self- 
energy factorji^iiii 



(f^)G^ (x,x') = -E^G^(x,x') 



(34) 



where E 



(imp) 



if) 



Ax +i/rjjt and the notations are used 



Q 



AJ =±Yv \ ^ 



U^fJ. 



dq Wiq + Xfj.) 



-oo 2vr I 2 - *l 



(35a) 



4) = 7c E IT [ W K* ~ x -) + W ^ + *k0 



(35b) 



Symbol V in (|35a|l stands for the integral principal value, 
the bar over the summation index in i|35bfl signifies that 
the summation is carried out over extended modes only. 
The conditional character of the term "impurity self en- 
ergy" with reference to S^ mp ' is related to the mere fact 
that this factor is actually determined by both the im- 
purity potential, whose correlator is proportional to the 
factor of Q, and the magnetic field, which renormalizcs 
the wavenumbers x^^ and also adjusts the number of 
extended modes, see next subsection. 

The action of the expressly non-local "magnetic" part 
of the operator (T^) is specified by the formula 



[ffi%](^ = E 



U^fJ. 



2iH 



S, 






Wll 



ly OX 



x ! dxi{G^{x,x x )) 



2iH 



C u 



' B 



WE 



l y dxi 



Ci>fJ. , 4 



H 2 
I'b 



G MM (a:i, x ) , 



(36) 



r 



wherefrom the "magnetic" self-energy, which is applica- 
ble under WS conditions, is immediately deduced 



S iB) 



-AH' 



E 












-f/x-jT" 



V *^ rfj, - *l 



>{u) 



>(u) 



(37) 



A. Mode content of the open quantum system 

Both the impurity and the magnetic self-energies are 
complex- valued quantities, whose real parts renormalize 
mode energies whereas the imaginary parts determine the 
uncertainty of energy levels. The requirement for mode 
energies to be positive defined specifies the number of ex- 
tended modes in the quantum system, which is normally 
referred to as the number of conducting channels, iV c . 
Computation of exact number of these modes, though 



clear in principle, is an intricate problem in general. For 
the system under consideration the number N c can be 
most easily found in the particular case of the magnetic 
field oriented lengthwise with respect to the current di- 
rection, i. e. for B || Ox. In this case mode energy 
renormalization due to the inter-mode magnetic scatter- 
ing, which is covered by the magnetic self-energy 1)37(1 . 
is small as compared with the intra-mode magnetic cor- 
rection present in the mode energy (|f7J) . Taking account 
of this fact, one can calculate the number of extended 
modes as 



N w 



N r 



E^V) 



where 



M z) 



hit 



k F H 



V1RI 



(38a) 



(38b) 



is the number of quantization levels in z-direction, whose 
energies lie beneath the Fermi energy, and 



N^\m) =int 



k F W 



k F H 



H 2 
Y1RI 



(38c) 



is the number of y-directional quantization levels perti- 
nent to the m-th level of z-quantization. Symbol int[. . .] 
in l|38b|) and (|38c|) denotes the integer part of the number 
enclosed in square brackets. 

The sum (|38a[) can be readily evaluated in the case 
where the number of extended modes relating to both 
of the transverse axes of the quantum waveguide is large 
as compared to unity. By replacing the sum with the 
integral one readily gets 

H 2 \ 
121%) ' 



N r 



47T 



1 



(39) 



wherefrom it is evident that application of the in-plane 
magnetic field can significantly reduce the number of ex- 
tended modes, even though inequality i|27[l holds true. 
This reduction is definitely the geometrical effect which 
is due to the curving of the electron orbits in the mag- 
netic field, and thus it can be only taken into considera- 
tion within the model of a finite- width quantum well that 
forms a 2D system. 

In Fig. |21 the numerical results for the number of 
effective conducting channels as a function of the in- 
verse magnetic field scaled as the Landau filling factor 



25 r 



15 - 



k F W/7r=10.5 


k F H/7r=5 
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k r H/7T=l 
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FIG. 2: The magnetic-field dependence of the number 
of effective conducting channels. The parameter v is the 
Landau filling factor, v — (fci^s) 2 - 

v = (kpls) 2 = kpR c are presented. The collapse of 
the number of current-carrying modes with a growth in 
the magnetic field is apparent, regardless of the quan- 
tum waveguide thickness H, the width W is assumed 
constant. The in-plane rotation of the magnetic field 
smoothly changes the presented picture because the real 
part of self-energy IJ37I) can at most reach the same (on 
the order of magnitude) value as the intra- mode magnetic 
addend in ljT7fl i22i 

In Fig. [31 the relation between the number of channels 
and the effective thickness of the quantum waveguide is 



presented, which actually demonstrates the dependence 
of N c on the depletion voltage adjusting the width of the 
near-surface potential well. In the extremely low mag- 
netic field (black curve) the number of channels increases 
nearly linear with growing H, in accordance with stan- 
dard geometrical consideration applicable to systems of 
waveguide configuration, and also with the conventional 
Ohm's law which is undoubtedly valid for bulk conduc- 
tors. With the growing magnetic field, the conventional 




FIG. 3: The number of conducting channels vs the 
width of the near-surface potential well at different val- 
ues of the Landau filling factor. The shaded area below 
the dashed curve covers the parameter region where the 
criterion (I27H is violated. 



geometric increase in the number of channels gets slower, 
gradually indicating the trend for lowering the number of 
conducting modes. This unusual dependence of the mode 
content of the electron waveguide on its effective thick- 
ness is due to non-monotonic dependence on H of the 
mode energy (|I7fl . 

Obviously, on a further increase of the magnetic field 
the tendency towards lowering the number of conduct- 
ing channels must be stabilized owing to terms oc lg 
in square brackets in r.h.s. of Eq. H3bj) . However, this 
can happen only in the domain of relatively strong mag- 
netic fields, where the WMS condition is violated and 
the approximate expression (|32|) for the inter-mode po- 
tential is no longer applicable. In such magnetic fields, 
the bulk quantum Hall effect is expected to come in the 
foreground, which is beyond the scope of this paper. 
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B. Dephasing of the mode states: the 
magnetic-field driven disorder 

Besides the impact on the number of extended quan- 
tum modes whose transverse energies are beneath the 
Fermi level, the in-plane magnetic field can significantly 
affect the coherent properties of the conducting chan- 
nels. This field controls the imaginary parts of both the 
impurity-governed self-energy l|35l) and the magnetic self- 
energy l|37|) . Both of these self-energies arise due to the 
inter-mode scattering. One should bear in mind, how- 
ever, that S^ m is basically determined by scattering 
from the impurity potential whereas the magnetic self- 
energy, Z 1 ^ , originates in the main from mode mixing 
due to the orbital effect of in-plane magnetic field. 

It is important to note that intermixing of channels 
which is controlled solely by the magnetic field cannot 
result in significant dephasing of mode states. By com- 
paring the imaginary part of self-energy l|37|) and the level 
width (|35b|) one can determine that the ratio of "purely 
magnetic" and "impurity-governed" dephasing rates is 
evaluated as 



QS, 



(B) 



-,(imp) 



H 



R c k F W 



H 



<l 



(40) 

This implies that under WMS condition the magnetic- 
field originated dephasing is negligible, whatever strength 
of the disorder. The conclusion is thus unavoidable that 
strong inter-mode mixing resulting from the magnetic 
field cannot give rise to significantly widening the mode 
levels unless there exists some random potential due to, 
say, impurities or the roughness of quantum well bound- 
aries, which can mediate the dephasing effect of the mag- 
netic field. The specific role of the magnetic field, as far 
as the mode entanglement is concerned, reduces to the 
change in collective parameters of the electron motion, 
such as the mode content of the confined system and 
the mode density of states, and in such an indirect way 
to modification of scattering parameters pertinent to ran- 
dom generators of inter-mode transitions (i. e., the impu- 
rity scattering cross-section, the polar pattern of electron 
reflection from rough boundaries, etc.). 

The influence of the magnetic field upon transport pa- 
rameters manifests itself directly through the mode de- 
phasing rate. Analytically, the estimate of this quantity 
can be most easily deduced from Eq. I)35bj) in the case 
where the number of quantization levels related to both 
of the transverse directions is large as compared to unity 
and the sum in Eq. (|35b|) can be replaced with the in- 
tegral. The dephasing rate for the particular mode fj, in 
this case reads 



1 



1 



» 



(B) 



» 



(0) 



12i?? 



(41) 



where 1/t^(0) = kpQ/iiT is the /n-th mode level width 
attributed to scattering due to the disorder potential 



only, with no external magnetic field^ The value of this 
zero-field level width equals exactly half the inverse mean 
free time calculated within the framework of classical ki- 
netic theory. Note that in the domain of weak magnetic 
fields corresponding to inequality (|27() the dephasing rate 
(|41|l decreases nearly quadratically in the magnetic field 
and has universal value, the same for each of the extended 
modest 

The result l|41|) . which is actually semiclassical, is of 
limited applicability. Upon varying the magnetic field the 
number of extended modes changes stepwise. Therefore 
the majority of physical quantities are bound to exhibit 
the oscillatory behaviour, which is closely related to well- 
known van Hove singularities in MDOS. In Fig. 0J the 
dephasing rates obtained numerically from Eq. (|35b|) for 
two specific modes of the electron waveguide are shown as 




5 v 



FIG. 4: The dephasing rate l|35b(l for two particular 
modes vs inverse magnetic field. 

functions of the inverse magnetic field. Square-root sin- 
gularities manifestly develop on both of the curves. One 
can also notice that scattering frequencies for different 
modes start to noticeably deviate from one another only 
in the range of relatively strong magnetic fields, where 
the number of extended modes assumes the value com- 
parable with unity. 

Besides the magnetic-field singularities depicted in 
Fig. QJ in Fig. the dephasing rate of the particular 
mode vs size parameters of the quantum waveguide is 
presented for two distinct values of |B|. Here, MDOS 
oscillations caused by abrupt changes in the number of 
conducting channels also make themselves very evident. 
We are led to conclude that by means of the orbital cou- 
pling to the electrons in a Q2D conducting system the 
in-plane magnetic field can take an effect which in some 
sense is analogous to that of electrostatic confinement 
potential. At the same time, in contrast to the magnetic- 
field-controlled singularities of the dephasing rate, which 
are depicted in Fig. 01 oscillations of truly geometrical ori- 
gin are noticeably more complicated. The distinction is 
caused by substantially different response of the effective 
mode energy (I17|l to the magnetic field, on the one hand, 
and to size parameters of the confined electron system, 
on the other. However, it should be noted that in both 
of the graphs, 4 and 5, the reduction of the dephasing by 
quenched disorder is clearly visible as the magnetic field 




1 1.5 2 2.5 3 3.5 ^/^ 

FIG. 5: The dephasing rate vs the quantum waveguide 
thickness at different strengths of in-piane magnetic fieid. 
The broken fraction of the lower curve falls into the range 
of parameters where WMS condition 1271 is violated. 



grows. This fact can serve as the indication of increas- 
ing coherence of electron transport in quench-disordered 
Q2D systems if they are subjected to external magnetic 
field. 



V. CONCLUSION 

In this study we have demonstrated that the observed 
giant positive magnetoresistance of 2D electron and hole 
systems subject to parallel magnetic field can be reason- 
ably explained in the framework of Fermi liquid theory 
being applied to structures created by confining poten- 
tial wells of finite rather than zero width. The magnetic 
field coupling to the carrier orbital motion which is due 
to finite thickness of guasi-two-dimensional layers, even 
though rather weak from semi-classical point of view, has 
been proven to influence quite essentially the collective 
electron spectrum. The reduction in the number of ex- 
tended modes with a growth of the magnetic field, as seen 
from Fig. |2 is very significant, continuing right up to zero 
in moderately strong fields, whereas individual electron 
trajectories in the plane normal to the magnetic field can 
go far beyond the effective thickness of the gated carrier 
system. The mode truncation effect of the in-plane mag- 
netic field is the more noticeable the larger is the aspect 
ratio of the confining potential well forming the electron 
waveguide. 

Since the number of extended modes, according to the 
Landauer theory, specifies the conductance of a bounded 
system, the results presented in Figs. and [3] can be di- 
rectly related to the experiment. In fact, they may be 
regarded as showing the conductance dependence on the 
corresponding parameters in the case of a perfect confin- 
ing potential well. As the perfect we mean a waveguide- 
type structure in which any mechanism of collective scat- 
tering of properly defined carrier modes does not exist. 
This actually implies that no scattering fields other than 



those involved in the unperturbed quasi-particle state 
formation in a particular system are taken into account. 
Specifically, the collective states pertinent to the problem 
considered in this study are specified by the confining 
potential profile. In the absence of the disorder poten- 
tial the collective electron motion should be regarded as 
ballistic, even though individual carriers do experience 
strong (specular) scattering at side boundaries of the po- 
tential well. 

If some random potential is involved, e. g., impurities 
or the roughness of quantum well boundaries, it should 
lead to stochastic rather than regular scattering of the 
primordial quasi-particles. It seems advantageous to sep- 
arate this type of scattering into two kinds, namely, intra- 
and inter-mode scattering. The former type of scatter- 
ing provides renormalization of transport parameters and 
also gives rise to Anderson localization of carrier states 
in the direction of current. The latter type, inelastic in 
form from the viewpoint of mode theory, leads to stochas- 
tic spreading of mode energy levels, or, in other words, 
to spatial dephasing of mode states. At first glance, it 
may appear that inter-mode scattering caused exclusively 
by the magnetic field is bound to produce the dephasing 
effect analogous to that introduced by the quenched dis- 
order. However, the estimate (|40f) is obviously contra- 
dicting to this expectation. According to the evaluation, 
in the absence of random potential, which ensures proba- 
bilistic property of mode energy levels, no imaginary part 
must be contained in the mode self-energy, in spite of 
substantial inter-mode mixing due to the magnetic field. 

Physically, this fact seems to be quite natural. Indeed, 
if one chooses to model lateral confinement of a Q2D car- 
rier system by the quadratic rather than the rectangular 
potential, eigen-functions of the transverse Hamiltonian 
could be obviously selected so as to completely avoid the 
mode coupling due to the magnetic field. The additional 
random potential, though static, would be in this case 
the only cause of the mode levels widening. At the same 
time, the quadratic confinement possesses the same sym- 
metry of the confined system as the rectangular well does. 
Therefore, it would be difficult to substantiate the drastic 
difference of the results obtained within the framework 
of these two models if one is guided by general consider- 
ations only. 

Fortunately, the result (|40l) reveals the lack (in the 
asymptotic sense) of the magnetic-field-originated de- 
phasing of the natural carrier spectrum. Certainly, the 
magnetic field does take part in the mode level spread- 
ing, yet mostly through the dependence on this field of 
the number of extended modes and of the mode den- 
sity of states. This parameters essentially determine 
the impurity-originated dephasing rate (|35b|l . which can 
thus be viewed as being produced by the magnetic-field- 
dependent disorder. The idea of the "magnetic-field- 
driven disorder" was previously suggested in Ref. |23, so 
the result l|35b|l can be viewed as substantiating the ra- 
tionality of such an interpretation. Clearly, in order to 
make a detailed comparison with experimental observa- 
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tions it is necessary to derive required formulas for the 
magnetocondactance. This work will be postponed for 
the next publication. 
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APPENDIX: DISORDER AVERAGING OF THE 
TRIAL GREEN FUNCTION 

After substitution of functions l|22J) into H20|l . the trial 
Green function inside magnetically biased interval C can 
be represented as a sum of four packets of spatial har- 
monics, viz. 



G (y) (x,x') =Qx(x,x')e i * { - x - x ' ) +g 2 (x,x')e- i * (x - x ' ) +g 3 (x,x')e i * (x+x " ) +g 4 (x,x')e-^ x+x ' ] 
Here, smooth envelope functions are given as 



Gi(x,x') = —A(x) 

z>c 

g 2 (x,x') = ^A(x) 

Zx 



e+ M^- e -^^y {x)e 

e _^l-e + r + (x)^^e 2 ^ 



7T+(x) 



7T_ (x) 



g 3 (x, x') = —A{x)J« L 

ZK 

gi{x.x') = -^-A{x)e^ L 

ZK 



7_(x') n+(x') ' 

B+ r~r + B_ r - r T^(x) 



e_ 



7T_ (x) 

n + (x) 
I 



n+(x) 
+ Q+T + (x) 



7T_(x') 



7T_(x) 



(A.l) 

(A.2a) 
(A.2b) 
(A.2c) 
(A.2d) 



where the notations are used 

A(x) = [i + r+(x)r_(x)e 2 ^ L ] 

r ± (x) =J±(x)/tt ± (x) , 
0± = 9[±{x-x')} . 



(A.3a) 
(A.3b) 
(A.3c) 



As regards the functions r±(x), their physical mean- 
ing is readily deduced from Eq. (|22l) . They represent 
reflection factors of spatial harmonics ±x incident at the 
point x onto the layers with end coordinates x and ±L/2, 
respectively. This factors meet the Riccati-type dynamic 
equations, 

±^T- = 2iv(x)r±(x) -c ± (z) + C£(z)r 2 ± (z) , (A.4) 

with boundary conditions stemming from l|25|) . 

Y±(±L/2) = %^ . (A.5) 

The averaging technique for the Junctionals of random 
fields (J23J was elaborated in Refs. H5H22H2I Here we 
only briefly indicate the main peculiarities of dealing with 
functionals of such a sort and present the result of the 
function l|A.l|) averaging. 

Having regard to correlation relations l|28|) it was 
p roven i5 J 22 i 23 f.]-^ binary correlation functions of the ef- 
fective random fields l|24|) under WIS conditions can be 



cast to the form 

(t]{x)t)(x')) = 

(C ± (a0Ci(O) 



■00 



Fi{x-x') 



Fi(x-x') 



-(V) 



(A.6a) 
(A.6b) 



where Lf and Lf, are the forward and the backward scat- 
tering lengths given in (|3l7|l for the particular mode fi. 
The function Fi(x) has the form 



™-£&-i$-*{ 1 - 



■21 



9(21 -\x\) 

(A.7) 

and plays the role of under-limiting (^-function when av- 
eraging smooth factors similar to the envelopes (|A.2(I . 
Before averaging the function (|A.1J| it makes sense to go 
over from functions 7±(x), 7r±(x) and C±( x ) to phase- 
renormalized functions 



l±(x) 


= 7-1- (x) exp 


,±L/2 

±i / dxi77(xi) 

J X 


(A.8a) 


7T±(x) 


— ir± (x) exp 


r±L/2 
Ti / dxir?(xi) 

J X 


, (A.8b) 


C±(s) 


= C± ( x ) ex P 


,±L/2 

±2i / dxi77(xi) 

Jx 


, (A.8c) 
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which enables one to remove the forward-scattering ran- 
dom field r)(x) from all dynamic equations and to sepa- 
rate it out in the form of exponential factors. In partic- 
ular, note then that correlation relation l|A.6b|l remains 
unchanged after renormalization ljA.8cf> . 

One can easily reveal that in view of short-range cor- 
relation of random functions i|24|) and due to the causal 
nature of functionals being averaged, the averaging of 
functionals with different sign indices in l|A.2|) can be 
done separately. By averaging the equation 



In view of smallness of the reflection coefficient Ji^ this 
allows one, when averaging l|A.l(l . to retain in l|A.2|) only 
the terms which do not contain factors r±(x) and 7±(x). 



± 



dT±(a;) 
da; 



-C ± (aO + C±(*)I±(aO (A.9) 



using the Furutsu-Novikov formula for gaussian random 
process^ we obtain 



In order to average the ratio tt±(x')/tt±(x), which is 
present in the principal terms of (|A.2fl . it is worthwhile 
to consider its Fourier transform over x' which, in view 
of the presence of 8-functions in l|A.2(l , takes the form 



(r ± (x)) = ^ B ) 



exp 



1 



Tx 



(A.10) 



±L/2 



& ± \x,q)=± I " dxi^^exp 

Jx K±(X) 



—iq{x — Xi)+in\x — Xi\±i I Ax2rj{x2) 



where forward-scattering random field rj(x) is already singled out. The averaging over this field yields 

\x-xi[ 



cxp 



±i I dx2rj(x 2 ) 

xx 



exp 



n \ 2L f 

and the function (|A.11|) . averaged beforehand over ry(x), is found to obey the equation 



T 



da; ~~ \2L, 



ixTiq) (<& (±) (x, q)) n ~ C± (x)T± (s)<* (±) (x, q)) 



(A.ll) 



(A.12) 



(A.13) 



which is to be solved along with Eq. (|A.9(1 . By averaging IJA.13JI over the effective back-scattering field £± (x) we arrive 
at the dynamic equations 



T 



d(<P^(x,g)) =1 



da- 



2 \L f + L b 



ix^f iq 



(*<*> (*,«)) 



with obvious "initial" conditions ($^(±£/2,g)) = 0. The solution to Eq. I|A.14[) has the form 

' L 



finally yielding 



\\T f + h '""'' 



1 — exp 



\\T S + T b ] " X{ l,J 



(g^x^'))^-^ + (g 2 (x,x'))e'^ x -^ 



2x 



exp 



1 ( 1 



2 \L 



1 

u 



2 TX 



\x — x 



(A.14) 



(A.15) 



(A.16) 



The envelopes (|A.2cj) and (|A.2dJl can be averaged in small in the parameter (|27|l and can thus be omitted, 
the same manner. Because both of them are proportional leaving the result (I A. 16|) as the main approximation for 
to reflection coefficient 7±(x), they prove to be relatively the impurity- averaged trial Green function. 
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